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We show that the Fourier transform of the Deeply Virtual Compton Scattering (DVCS) am- 



^ ■ plitude with respect to the skewness variable C at fixed invariant momentum transfer squared t 

provides a unique way to visualize the structure of the target hadron in the boost-invariant longi- 
tudinal coordinate space. The results are analogous to the diffractive scattering of a wave in optics. 
As a specific example, we utilize the quantum fiuctuations of a fermion state at one loop in QED 
to obtain the behavior of the DVCS amplitude for electron-photon scattering. We then simulate 
the wavefunctions for a hadron by differentiating the above LFWFs with respect to M^ and study 
the corresponding DVCS amplitudes in light-front longitudinal space. In both cases we observe 
that the diffractive patterns in the longitudinal variable conjugate to ( sharpen and the positions 
of the first minima move in with increasing momentum transfer. For fixed t, higher minima appear 
at positions which are integral multiples of the lowest minimum. Both these observations strongly 
support the analogy with diffraction in optics. 



Introduction 

Deeply Virtual Compton Scattering (DVCS) 7*(g) + p{P) -^ 7(?') + p(-P') provides a 
remarkable tool for studying the fundamental structure of the proton at the amplitude level. 
We define the momentum transfer A = P~P' and the invariant momentum transfer squared 
t = A'^A^. When the incoming photon is highly virtual Q^ = —q^ » Ag^^ , the underly- 
ing scattering process measures Compton scattering on bound quarks, convoluted with the 
fundamental microscopic wavefunctions of the initial- and final-state proton. In addition, 
the initial-state photon can scatter on virtual qq pairs in the target which are then annihi- 
lated by the final-state photon, thus probing the particle-number quantum fluctuations of 
the hadron wavefunction required for Lorentz invariance. Measurements of the DVCS cross 
sections with specific proton and photon polarizations can provide comprehensive probes of 
the spin as well as spatial structure of the proton at the most fundamental level of QCD. 

The theoretical analysis of DVCS is particularly clear and compelling when one utilizes 
light-front quantization at fixed t = y^. (We use the standard LF coordinates P^ = 
po _|_ p'^^y^ = y^ ± y^. Since the proton is on-shell, P^P~ — P± = M^). If we neglect 
radiative corrections to the struck quark propagator (i.e., set the Wilson line to 1), then 
the required DVCS quark matrix elements can be computed from the overlap of the boost- 
invariant light- front Fock state wavefunctions (LFWFs) of the target hadron [l|, |2]. The 
longitudinal momentum transfer to the target hadron is given by the "skewness" variable 
( = ^— . Since the incoming photon is space-like (g^ < 0) and the final photon is on-shell 
((g')^ = 0), the skewness is never zero in a physical experiment. The DVCS process involves 
off-forward hadronic matrix elements of light-front bilocal currents. Accordingly, in different 
kinematical regions, there is a diagonal parton-number conserving n ^ n overlap and an 
off-diagonal n + 1 ^ n — 1 overlap where the parton number is decreased by two. Thus, 
given the LFWFs one then obtains a complete specification of all of the Generalized Parton 
Distributions (GPDs) measurable in DVCS, including theirphase structure. The sum rules 
of DVCS, such as Ji's sum rule for angular momentum [j] and the integral relations to 
electromagnetic and gravitational form factors are all explicitly satisfied in the light-front 
(LF) formahsm [J,!^. 

In this paper we shall show how one can use measurements of the dependence of the DVCS 
amplitude on the skewness variable ( to obtain a novel optical image of a hadron target, in 



analogy to the way in which one scatters optical waves to form a diffraction pattern. In this 
context, it is useful to introduce a coordinate b conjugate to the momentum transfer A such 
that 6 • A = I&+A- + |6"A+ - b± ■ A±. Note that |&"A+ = |&"P+C = o^C where we have 
defined the boost invariant variable a which is an 'impact parameter' in the longitudinal 
coordinate space. The Fourier transform of the DVCS amplitude with respect to ( allows 
one to determine the longitudinal structure of the target hadron in terms of the variable a. 

Burkardt j^ |5[ has studied the off-forward parton distribution function at zero longitudi- 
nal momentum transfer and fixed longitudinal momentum fraction x in the impact parameter 
(6_l) space. This relativistic impact representation on the light-front was introduced earlier 
by Soper 'Jo] in the context of the Fourier Transform of the elastic form factor. We study the 
DVCS amplitude, which involves integration over x, in the a space at fixed four-momentum 
transfer —t. Note that experimentally DVCS amplitudes are measured as a function of C 
and —t. Thus, our work is suited for the direct analysis of experimental data and is com- 
plementary to the work of Burkardt and Soper. If one combines the longitudinal transform 
(at fixed A_l) with the Fourier transform (FT) of the DVCS amplitude with respect to A_l 
one can obtain a complete three-dimensional description of hadron optics at fixed LF time. 

Recently, a 3D picture of the proton has been proposed in f^ in a different approach, 
in terms of a Wigner distribution for the relativistic quarks and gluons inside the proton. 
A major difference of this from our work is that the Wigner distributions are defined in 
the rest frame of the proton. Integrating out k~ one gets the reduced Wigner distributions 
which are not observable quantities in the quantum domain. Further integration over k^ 
relates them to the FT of GPDs H{x, C,, t) and E{x, ^, t), where ^ = q^ /2Eq and x is a special 
combination of off shell energy and momentum along z. On the other hand, we are taking 
the FT of the experimentally measurable DVCS amplitudes directly and not the GPDs. 

In principle, the LFWFs of hadrons in QCD can be computed using a nonperturbative 
method such as Discretized Light Cone Quantization (DLCQ) where the LF Hamiltonian 



is diagonalized on a free Fock basis |S|. This has been accomplished for simpler confining 
quantum field theories such as QCD{1 + 1) 9]. Models for the LFWFs of hadrons in 
(3+1) dimensions displaying confinement at large distances and conformal symmetry at short 



distances have been obtained using the AdS/CFT method [10|. The 



Q. 



space structure of topological objects has been studied in DLCQ 
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In order to illustrate our general framework, we will present here an explicit calculation 



of the a transform of virtual Compton scattering on the quantum fluctuations of a lepton in 
QED at one-loop order [12] , the same system which gives the Schwinger anomalous moment 
ajl-n . This model has the advantage that it is Lorentz invariant, and thus it has the correct 
relationship between the diagonal n ^ n and the off-diagonal n — 1 ^ n + 1 Fock state 
contributions to the DVCS amplitude. One can generalize this analysis by assigning a mass 
M to the external electrons and a different mass m to the internal electron lines and a mass 
A to the internal photon lines with M < m + \ for stability. In effect, we shall represent 
a spin-| system as a composite of a spin-| fermion and a spin-1 vector boson [13|, |lj, |l5[- 
We also will present numerical results for a composite hadron by taking a derivative of the 
LFWFs with respect to the hadron's mass M^. This simulates the behavior of a bound-state 
hadron by improving the fall-off at the end points of the longitudinal momentum fraction 
X. The summary of our main results will be given in this letter. A more detailed analysis 
will be given in a forthcoming article |lq . 

DVCS in the LF Formalism 

The kinematics of the DVCS process has been given in detail in [jj, |^. One can work 
in a frame where the momenta of the initial and final proton has a A -^ —A symmetry 
121 . However, in this frame, the kinematics in terms of the parton momenta becomes more 

n 

complicated. Here, we choose the frame of Ref. [1||. 

The virtual Compton amplitude M'^'^(gj_, A_l, (^), i.e., the transition matrix element of 
the process 7*(g) + p{P) — » 7(?') +p(P'), can be defined from the light-cone time-ordered 
product of currents 

M^'^(g-1,A^,C) = ijdSe-''^-y{P'\T.ny)rmP) , (1) 

where the Lorentz indices /i and v denote the polarizations of the initial and final photons 
respectively. In the limit Q"^ ^ oo aX fixed C, and t the Compton amplitude is thus given by 

M"(g-1,A^,C) = e^ef M'^'^(g-1,A^,C) = - 
X I t'\z, C) U{P') [h{z, C, t) 7+ + E{z, C, t) 

where P = |(-P' + P) and we take a frame in which g+ 
compared to the masses and —t, we have, 
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up to corrections in 1/Q'^. Thus C, plays the role of the Bjorken variable in deeply virtual 
Compton scattering. For a fixed value of — t, the allowed range of C, is given by 



< C < 
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For simplicity we only consider one quark with fiavor q and electric charge e^. We here 
consider the contribution of only the spin-independent GPDs H and E. Throughout our 
analysis we will use the "handbag" approximation where corrections to the hard quark 
propagator are neglected. 

For circularly polarized initial and final photons (/, J are t or |)) contributions only 
come from 
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For a longitudinally polarized initial photon, the Compton amplitude is of order 1/Q and 
thus vanishes in the limit Q^ —>■ oo. At order 1/Q there are several corrections to the simple 
structure in Eq. (0). We do not consider them here. 

The generalized parton distributions H, E are defined through matrix elements of the 
bilinear vector and axial vector currents on the light-cone: 

dy- 
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The off-forward matrix elements given by Eq. (jU)) can be expressed in terms of overlaps of 
LFWFs of the state [l|, |2| • For this, we take the state to be an electron in QED at one loop 
and consider the LFWFs for this system. 

The light-front Fock state wavefunctions corresponding to the quantum fiuctuations of a 
physical electron can be systematically evaluated in QED perturbation theory. The state is 
expanded in Fock space and there are contributions from | e~7) and | e~e~e~^), in addition 
to renormalizing the one-electron state. The two-particle state is expanded as, 

dxd^A; 
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where the two-particle states \Sf,s^; x,k±) are normahzed as in [1|. The variables Sf and 
s^ denote the projection of the spins of the constituent fermion and boson along the quan- 
tization axis, and the variables x and k± refer to the momentum of the fermion. The light 
cone momentum fractions Xi = ^ satisfy < Xi < 1, J2i^i = 1- We employ the light-cone 
gauge A'^ = 0, so that the gauge boson polarizations are physical. The three-particle state 
has a similar expansion. Both the two- and three-particle Fock state components are given 
in [l| . We list here the two-particle wavefunctions for the spin- up electron l|, |lj, |l7 1 
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Similarly, the wavefunction for an electron with negative helicity can also be obtained. 

In the domain ^ < 2; < 1, there are diagonal 2^2 overlap contributions to Eq. ©, both 
hehcity flip, Ff_ (A' ^ A) and helicity non-flip, Ff+ (A' = A) [1]. The GPDs H(^2^2){zX,t) 
and -E(2-^2)(-2; C)^) ^^e zero in the domain ( — 1 < z < 0, which corresponds to emission 
and reabsorption of an e"*" from a physical electron. Contributions to if(„^„)(z, C, t) and 
E{n^n){z,Cjt) ill that domain only appear beyond one-loop level. This is because in the 
DVCS amplitude we have integrations over z, y~, and x. When integration over y~ is 
performed, the fermion part of the bilocal current yields a factor 6{z — x) and the anti- 
fermion part of the bilocal current yields a factor 6{z + x). The latter contribution is absent 
in the one loop DVCS amplitude of a electron target, which we consider in the present work. 

The matrix elements F^^ and F^^_ are calculated using the two-particle LFWFs given 
in Eq. (jHl). The contributions in the domain, < z < (, namely, F^^ and F^\_ come 
from overlaps of three-particle and one-particle LFWFs [l|. These are calculated using the 
three-particle wavefunction. Explicit expressions of all the above matrix elements will be 
given in jl6|. 



We calculate the DVCS amplitude given by Eq. (0) using the off-forward matrix elements 
calculated above. In order to regulate the ultraviolet divergences, we use a cutoff A on the 
transverse momentum fc-*-. The real and imaginary parts are calculated separately using the 



prescription 

r dx ^—-F{x, = P C dx^F{x, C) - inFiC, (). (10) 

Jo X — (, + te JO X — C 

Here P denotes the principal value defined as 

P /' dx^Fix, = lim [ /^"' -^F{x, + f' -^F{x, ()] (11) 

Jo X — Q e^OLJo X — Q Jc.+eX — C, ^ 

where 



F(x,C) = i^^-^(x,C,A^), forO<a;<C 
= i^f(x,C,A±), forC<a;<l 

with ij = ++ for helicity non-fiip and ij = H — for helicity fiip amplitudes. Since the off- 
forward matrix elements are continuous at x = C, -^(C)C) = Pf^i.^ = (,(^^±) = Pf^i.^ = 
C,C;A_l). Note that the principal value prescription cannot be used at x = 0. We take a 
small cutoff at this point for the numerical calculation. The off-forward matrix elements F^^ 
(which contribute in the kinematical region < x < () vanish as x — > 0, as a result there 
is no logarithmic divergence at this point for nonzero (. But, we need to be careful here as 
when we consider the Fourier transform in a space, ( can go to zero and divergences from 
small X can occur from F^^ which is finite and nonzero at x, C — ^ 0. 

The imaginary part of the amplitude when the electron helicity is not fiipped is then 
given by 

Im[M++](C, Ax) = ne'Pllix = C,C, A^). (12) 

A similar expression holds in the case when the electron helicity is fiipped (lm[M+_](C, A_|_)) 
in which F_li are replaced by F+^. The helicity-fiip DVCS amplitude is proportional to 
(Ai —iA2) [1&\. Without any loss of generality, the plots for these amplitudes are presented 
with A2 = 0. The imaginary part receives contributions at x = C. The other regions 
of X contribute to the real part. It is to be emphasized that we are using the handbag 
approximation of the DVCS amplitude. Contributions from the Wilson lines are in general 
not zero, and they can give rise to new phase structures as seen in single-spin asymmetries 

0. 

The real part of the DVCS amplitude in our model is given by 

C-^i .-, , . , r 1 

- + 



Re [M++] (C,Ax) = -e^^ '' dx Ff+(x,C, A 



X X — ( 



-e' f~' dx Pllix, C, Ax) [ - + ^ ] . (13) 
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A similar expression holds for the helicity flip DVCS amplitude. The cutoff dependence at 
X = C in the principal value prescription gets canceled explicitly and, as a result, the DVCS 
amplitude is independent of the cutoff ei . 

Calculation of the a Fourier Transform 

In order to obtain the DVCS amplitude in longitudinal coordinate space, we take a Fourier 
transform in (^ as, 

A+_(a,t) = i-£"''c?C e^'^^ M+_(C, Ax), (14) 

where a = ^P^b~ is the (boost invariant) longitudinal distance on the light-cone and the 
FTs are performed at a flxed invariant momentum transfer squared —t. We have imposed 
cutoffs at e = ci = 69/2 = 0.001 for the numerical calculation. 

n 

A detailed discussion of the cutoff scheme will be given in [16j . 

All Fourier transforms (FT) have been performed by numerically calculating the Fourier 
sine and cosine transforms and then calculating the resultant by squaring them, adding and 
taking the square root, thereby yielding the Fourier Spectrum (FS). In Fig. 1, we have shown 
the FS of the imaginary part of the DVCS amplitude for M = 0.51 MeV, m = 0.5 MeV and 
A = 0.02 MeV. (a) is the helicity non-flip and (b) is the helicity flip part of the amplitude. 
We have divided the amplitude by the normalization constant e^/(167r'^) and have taken 
A = Q = 10 MeV. The helicity non-flip amplitude depends on the scale A logarithmically 
and the scale dependence is suppressed in the helicity-flip part. As seen in Fig. 1(a), the FS 
of the imaginary part of the helicity non-flip amplitude displays a diffraction pattern in a. 
The peak initially increases with increasing —t, but then decreases as —t increases further. 
The latter behavior warrants further study. In contrast, we see in Fig. 1(b) that there is 
no diffraction pattern in the FS of the imaginary part of the helicity flip amplitude. This is 
due to different behavior with respect to ( of the respective amplitudes [16| . Further, note 
that in this case the peak monotonically increases with increasing —t. One reason for this 
may be the presence of the extra factor of A^ in the helicity flip amplitude compared to the 
helicity non-flip amplitude. 

In Fig. 2 we have plotted the Fourier Spectrum of the real part of the DVCS amplitude 
vs. cr for M = 0.51 MeV, m = 0.5 MeV and A = 0.02 MeV. For the same | t |, the behavior 
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FIG. 1: Fourier spectrum of the imaginary part of the DVCS amphtude of an electron vs. a for 
M = 0.51 MeV, m = 0.5 MeV, A = 0.02 MeV, (a) when the electron hehcity is not flipped; (b) 
when the helicity is flipped. The parameter t is in MeV . 
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FIG. 2: Fourier spectrum of the real part of the DVCS amplitude of an electron vs. a for M = 0.51 
MeV, m = 0.5 MeV, A = 0.02 MeV, (a) when the electron helicity is not flipped; (b) when the 
helicity is flipped. The parameter t is in MeV . 

is independent of Q when | t |< m'^. For each Q, the peak at cr = is sharper and higher 
as I t I increases, the number of minima within the same a range also increases. As seen in 
Fig. 2, the FS of both the helicity flip and helicity non-flip parts show a diffraction pattern. 



The DVCS amplitude for an electron-like state at one loop has potential singularities at 
X = 1. As mentioned above, we have used cutoffs at a; = 0, 1. The cutoff at s = is imposed 
for the numerical integration. In the 2— and 3-body LFWFs, the bound-state mass squared 
M^ appears in the denominator. Differentiation of the LFWFs with respect to M^ increases 
the fall-off of the wavefunctions near the end points x = 0, 1 and mimics the hadronic 
wavef unctions. In this way, the cutoff dependency is removed. Differentiating once with 
respect to M^ simulates a meson-like wavefunction and another differentiation simulates a 
proton wavefunction. Convolution of these wavefunctions in the same way as we have done 
for the dressed electron wavefunctions will simulate the corresponding DVCS amplitudes 
for bound state hadrons. One has to note that differentiation of the single particle wave 
function yields zero and thus there is no 3 — 1 overlap contribution to the DVCS amplitude 
in this hadron model. It is to be noted that in recent holographic models from AdS/CFT 
as well jlO|, only valence LFWFs are constructed. 

The equivalent but easier way is to differentiate the DVCS amplitude with respect to the 
initial and final state masses. Here we calculate the quantity Mp-Qj^Mj-^^AijiyMi^Mp) 
where M/, Mp are the initial and final bound state masses. For numerical computation, we 
use the discrete (in the sense that the denominator is small and finite but not limiting to 
zero) version of the differentiation 

where M^ = iMi+Mil and 6M^ = (M^ - M^). We have taken Mn, Mpi = 150 + 1, 
Mi2, Mf2 = 150-1 MeV and fixed parameters M = 150 and m = A = 300 MeV. In Figs. 3 
and 4 we have shown the DVCS amplitude of the simulated hadron model, both as a function 
of ( and after taking the FT in (. In Fig. 4 (c), we have plotted the structure function 
F2{x) in this model. The wave function is normalized tol. Recall that the 7*p —>■ jp 
DVCS amplitude has both real [19| and imaginary parts 20]. The imaginary part requires 
a non-vanishing LFWF at x' = jr^ = 0. If we consider a dressed electron, the imaginary 
part from the pole at x = ( survives because of the numerator -^^ factor in the electron's 
LFWF. This numerator behavior reflects the spin-1 nature of the constituent boson. The 
X — ( -^ singularity is shielded when we differentiate the final state LFWFs with respect 
to M^ and, as a result, the imaginary part of the amplitude vanishes in this model. We thus 
have constructed a model where the DVCS amplitude is purely real. It is interesting that the 
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FIG. 3: Real part of the DVCS amplitude for the simulated meson-like bound state. The param- 
eters are M = 150, m = X = 300 MeV. (a) Helicity flip amplitude vs. (, (b) Fourier spectrum of 
the same vs. a. The parameter t is in MeV . 

forward virtual Compton amplitude 7*p — > 7*p (whose imaginary part gives the structure 
function) does not have this property. The pole at a; = C is not shielded since the initial 
and final n = 2 LFWFs are functions of x. If instead we consider the differentiation with 
respect to the internal fermion mass m^ rather than the bound state mass M^, although it 
does not improve the wavefunction behavior at the endpoint x = 0, we can generate a model 
with both real and imaginary parts of the DVCS amplitudes. It is worthwhile to point out 
that in general the LFWFs for a hadron may be non- vanishing at the end points |2l|, and 
recent measurements of single spin asymmetries suggest that the GPDs are non-vanishing 
at x = C [22]. A more realistic estimate would require non- valence Fock states J2j|. 

From the plots, we propose an optics analog of the behavior of the Fourier Spectrum 
of the DVCS amplitude. In fact, the similarity of p araxial optics and quantum fields on 



the light cone was first explored long ago in 



2i 



25[. In the case of DVCS, the final-state 



proton wavefunction is modified relative to the initial state proton wavefunction because 
of the momentum transferred to the quark in the hard Compton scattering. The quark 
momentum undergoes changes in both longitudinal {(P^) and transverse (A^) directions. 
We remind the reader that, to keep close contact with experimental analysis, we have kept —t 
fixed while performing the Fourier transform over the ( variable. Note that the integrals over 
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FIG. 4: Real part of the DVCS amplitude for the simulated meson-like bound state. The param- 
eters are M = 150, m = X = 300 MeV. (a) Helicity non-flip amplitude vs. C, (b) Fourier spectrum 
of the same vs. cr, (c) Structure function vs. x. The parameter t is in MeV . 

X and ( are of finite range. More importantly the upper limit of ( integral is (max which in 
turn is determined by the value of —t. The finiteness of slit width is a necessary condition for 
the occurrence of diffraction pattern in optics. Thus when integration is performed over the 
range from to (max, this finite range acts as a slit of finite width and provides a necessary 
condition for the occurrence of diffraction pattern in the Fourier transform of the DVCS 
amplitude. When a diffraction pattern is produced, in analogy with single slit diffraction, 
we expect the position of the first minimum to be inversely proportional to (max- Since 
(max increases with —t, we expect the position of the first minimum to move to a smaller 
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value of 0", in analogy with optical diffraction. In the case of the Fourier Spectrum of DVCS 
on the quantum fluctuations of a lepton target in QED, and also in the corresponding 
hadronic model, one sees that the diffractive patterns in a sharpen and the positions of 
the first minima typically move in with increasing momentum transfer. Thus the invariant 
longitudinal size of the parton distribution becomes longer and the shape of the conjugate 
light-cone momentum distribution becomes narrower with increasing | t |. Regarding the 
diffraction patterns observed in the Fourier Spectrum of the DVCS amplitude, we further 
note that for fixed — t, higher minima appear at positions which are integral multiples of the 
lowest minimum. This further supports the analogy with diffraction in optics. 

We can study the diffraction pattern in cr as a function of t or A^ in order to register the 
effect of a change in transverse momentum resulting from the Compton scattering. If one 
Fourier transforms in C at fixed Ai and then Fourier transforms the change in transverse mo- 
■nentu. A. to „„pact space .. EtS. the„ o.e would have the analog of a th.ee-di.en.onal 
scattering center. In this sense, scattering photons in DVCS provides the complete Lorentz- 
invariant light front coordinate space structure of a hadron. 
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